Mass imbalance effect in resonant Bose-Fermi mixtures 
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We consider a homogeneous Bose-Fermi mixture, with the boson-fermion interaction tuned by a 
Fano-Feshbach resonance, in the presence of mass and density imbalance between the two species. 
We first study the finite temperature phase diagram for the specific case of the mass-imbalanced 
mixture 87 Rb — 40 K for different values of the density imbalance. We then analyse the quantum phase 
transition associated with the disappearance at zero temperature of the boson condensate above a 
critical boson-fermion coupling. We find a pronounced dependence of the critical coupling on the 
mass ratio and a weak dependence on the density imbalance. For a vanishingly small boson density, 
we derive the asymptotic expressions for the critical coupling in the limits of small and large mass 
ratios. These expressions are relevant also for the polaron-molecule transition in a Fermi mixture 
at small and large mass ratios. The analysis of the momentum distribution functions at sufficiently 
large density imbalances shows an interesting effect in the bosonic momentum distribution due to 
the simultaneous presence of composite fermions and unpaired fermions. 

PACS numbers: 03.75.Ss,03.75.Hh,32.30.Bv,74.20.-z 



I. INTRODUCTION 

An impressive series of experiments has been realized 
with ultracold atomic gases over the last years, reproduc- 
ing systems, or physical situations, relevant to several ar- 
eas of physics. The use of Fano-Feshbach resonances to 
control the interaction strength between particles, in par- 
ticular, has been the cornerstone of many of these recent 
experimental achievements. 

Besides providing analog quantum models and simula- 
tors for systems of interest in different fields of physics, 
ultracold gases offer also the possibility to construct novel 
many-body systems with no corresponding counterparts 
in other domains of physics. Resonant Bose-Fermi mix- 
tures constitute an interesting example in this respect, 
and have been the object of active theoretical and exper- 
imental investigation recently. 

Nonresonant Bose-Fermi mixtures were initially stud- 
ied theoretically and experimentally in [Tj-Q and re- 
spectively. Bose-Fermi mixtures in the presence of a nar- 
row Fano-Feshbach resonance were then considered in 0- 
Hfl |. For a narrow resonance one has to take into ac- 
count explicitly the molecular state forming in the closed 
channel, leading to a Hamiltonian which includes three 
different species (bosons, fermions, and molecules) from 
the outset. Calculations with this Hamiltonian have been 
done in the above works by using approximations strictly 
valid only in the limit of a very narrow resonance. This 
limit, however, does not seem to be applicable to most of 
current experiments on Bose-Fermi mixtures [l7l - [23j . 

A broad Fano-Feshbach resonance is characterized by 
the smallness of the effective range parameter ro of the 
boson-fermion scattering amplitude with respect to both 
the average interparticle distance and the boson-fermion 
scattering length a [24(. Under these conditions, the sys- 
tem can be described by a Hamiltonian made just by 
bosons and fermions mutually interacting via an attrac- 



tive contact potential. 

Initial works studying Bose-Fermi mixtures in the pres- 
ence of a broad resonance focused on lattice models [25r - 
[28j . or considered separable interactions, as inspired by 
nuclear-physics models [29(. The continuum case in the 
presence of an attractive contact potential was first tack- 
led in Ref. [30| . which concentrated on the thermody- 
namic properties of the condensed phase of a Bose-Fermi 
mixture at zero temperature. 

A first study of the competition between Bose-Fermi 
pairing and boson condensation in a broadly resonant 
Bose-Fermi mixture has been presented in a previous 
work by us [3l| . By using many-body diagrammatic the- 
ory, with a T-matrix approximation for the bosonic and 
fcrmionic self-energies, we were able to show that, for in- 
creasing Bose-Fermi attraction, the boson-fermion pair- 
ing correlations progressively reduce the boson condensa- 
tion temperature and make it eventually vanish at a crit- 
ical coupling above which the condensate is completely 
depleted (thus revealing the presence of a quantum phase 
transition at zero temperature). In the present paper 
we extend the work of Ref. [31] by analyzing the effect 
of a mass imbalance between the bosonic and fcrmionic 
species. We present results both at finite temperature 
(for a specific mass ratio) and at zero temperature (for 
several mass ratios). The effect of mass imbalance in 
a broadly resonant Bose-Fermi mixture has been studied 
recently also in Ref. (32j , within a path-integral approach 
which is complementary to our approach and that was 
limited to zero temperature only. A comparison between 
our results and those of Ref. [32| will be discussed later 
on in our paper. 

The paper is organized as follows. In section |TT] we 
present the theoretical formalism and the main equa- 
tions describing the physical system of interest in our 
paper. Section IIIII reports the numerical results at fi- 
nite temperature for a 87 Rb — 40 K mixture. The phase 



diagram for different values of the density imbalance is 
presented, together with the curves for the boson and 
fermion chemical potentials at the critical temperature. 
Section ITVl treats the zero temperature limit. We report 
the results for the dependence of the critical coupling on 
the mass ratio and density imbalance, the correspond- 
ing chemical potentials, and the momentum distribution 
functions. Section [V] presents finally our conclusions. 

II. FORMALISM 

We consider a homogeneous mixture, composed by 
single-component fermions and bosons. The boson- 
fermion interaction is assumed to be tuned by a broad 
Fano-Feshbach resonance, as in most of current experi- 
ments on Bose- Fermi mixtures [13-111] • Under this con- 
dition, the boson-fermion interaction can be adequately 
described by an attractive point-contact potential. We 
consider then the the following (grand-canonical) Hamil- 
tonian: 

+ v JdriPl(T)4(r)Mr)Mr)- (1) 

Here ^J(r), creates a particle of mass m s and chemi- 
cal potential fi s at spatial position r, where s=B,F indi- 
cates the boson and fermion atomic species, respectively, 
while vo is the bare strength of the contact interaction 
(we set h — fee = 1 throughout this paper). As for two- 
component Fermi gases [3J. the ultraviolet divergences 
associated with the contact interaction in (JTJ) are elimi- 
nated by expressing the bare interaction vq in terms of 
the boson-fermion scattering length a: 

1 m r f G?k 2m r 

m = 2^ ~ J (2tt) 3 k 2 ' U>) 

where m r — m^m-p/{m^, + mp) is the reduced mass of 
the boson-fermion system. 

We have not considered in the Hamiltonian (|T|) an ex- 
plicit boson-boson interaction term. In the physical sys- 
tems relevant to experiments, the corresponding nonres- 
onant scattering length is normally small and positive. 
In the homogeneous and normal system we are going to 
consider, it yields then only a mean- field shift of the bo- 
son chemical potential. Note however that, according 
to the variational analysis of Ref. [Hj], the assumption 
of having a homogeneous system actually requires the 
boson-boson scattering length to exceed a certain value, 
in order to guarantee the mechanical stability of the Bosc- 
Fermi mixture against collapse. We thus implicitly as- 
sume to work in this stable regime. Note finally that 
a Fermi-Fermi s— wave scattering length is excluded by 
Pauli principle. 

A natural length scale for the system is provided by the 
average interparticle distance n -1 / 3 (where n = riB + n F 
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FIG. 1: (Color online) T-matrix diagrams for the fermionic 
and bosonic self-energies in the normal phase. Full lines repre- 
sent bare bosonic (BB) and fermionic (FF) Green's functions. 
Broken lines represent bare boson-fermions interactions vo. 

is the total particle-number density, 7Jb and 7ip being 
the individual boson and fermion particle-number den- 
sity, respectively). We thus introduce a fictitious Fermi 
momentum of the system fcp = (3n 2 n) 1 ^ 3 (as for an 
equivalent two-component Fermi gas with a density equal 
to the total density of the system), and use the dimen- 
sionlcss coupling parameter g = (fcpa) -1 to describe the 
strength of the interaction. In the weak-coupling limit, 
where the scattering length a is small and negative and 
g <C —1, the two components behave essentially as Bose 
and Fermi ideal gases. In the opposite, strong-coupling, 
limit where a is small and positive and g » 1, the sys- 
tem is expected to be be described in terms of composite 
fermions, i.e. boson-fermion pairs with a binding energy 
€q = l/(2m r a 2 ), and excess fermions. [We restrict our 
analysis to mixtures where the number of bosons never 
exceeds the number of fermions. This is because only in 
this case a quantum phase transition associated with the 
disappearance of the condensate is possible [31| . In addi- 
tion, mixtures with riB > wp are expected to be severely 
affected by three-atom losses.] 

The physical behavior of a Bose-Fermi mixture across 
all the resonance is captured by the T-matrix set of di- 
agrams for the boson and fermion self-energies repre- 
sented in FigfTJ As shown in our previous work [3l| , this 
set of diagrams is able to recover the expected physical 
behavior in both the weak- and strong-coupling limits, 
thus providing a meaningful theoretical framework for 
the whole resonance. The corresponding equations for 
the bosonic and fermionic self-energies Sb and Ep and 
many-body T— matrix T read: 

(^E G F( p -<?)r(P) (3) 

SfW = T f^J2 G °B(P~k)T(P), (4) 

r <-- n ->~{5; + /<& 
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/Kf(P-p)] + 6[€b(p)] 2m r 



L Cf(p-p) + Cb(p)-^ 



(5) 



Here q = (q.,uj v ), fe = (k,w„), P = (P,Q m ), where 
cj„ = 2-kvT and u;„ = (2n + l)7rT, O m = (2to + 1)ttT 
are bosonic and fcrmionic Matsubara frequencies, respec- 
tively, (is,n,m being integer numbers), while f(x) and 
b(x) are the Fermi and Bose distribution functions at 
temperature T, and £ s (p) = P 2 /(2m s ) — fJ, s - 

The self-energies (|3]) and ((4]) determine the dressed 
Green's functions G s via the Dyson's equation G.T 1 = 
G° _1 — E s (with the bare Green's functions given by 
G° (q) = [iu v - ^(q)]- 1 and G° F (k) = [iuj n - ^(k)]- 1 ). 
The dressed Green's functions G s allow to calculate 
the boson and fermion momentum distribution functions 
through the equations: 

n B (q) = -r^G B (q,^)e l ^ 0+ (6) 
n F (k) = T^G F (k, W „) e ^" 0+ , (7) 



for the bosonic and fermionic self-energies, where q = 
(q, ub), k = (k, wp), P = (P,f2) while the frequencies 
wb, wp and fi are now continuous variables. 

Similarly, the equations ([5]) and ([7]) for the momentum 
distribution functions are changed to 



"B(q) = - j^-G B ( q ,UJ B )e^ 



np(k) 



~2V 



(13) 
(14) 



from which the number densities can be calculated as 
before. 

A closed-form expression can be finally derived for the 
many-body T-matrix r(P) at zero temperature, when 
the Fermi and Bose distribution functions appearing in 
Eq. ([5]) are replaced by step functions. We report in par- 
ticular the expression for T(P) when uf > and /ib < 0, 
which is in practice the only one relevant to our calcula- 
tions at zero temperature. 

We have 



which in turn determine the boson and fermion number 
densities: 



7lF 



1 1 



n B (q) 



(2tt) 3 
dk 

(2^ nF(k) - 



(8) 
(9) 



At fixed densities and temperature, the two coupled 
equations © and ([§]) fully determine the boson and 
fermion chemical potentials, and therefore the thermody- 
namic properties of the Bose- Fermi mixture in the normal 
phase. 

Coming from the normal phase, the condensation of 
bosons starts when the condition 



M B-S B (g = Q) = 



is first met. Eq. (|T0|) corresponds to the requirement of 
a divergent occupancy of the boson momentum distribu- 
tion at zero momentum: lim q _ ! .o ^B(q) = oo. 35] 



A. The zero-temperature limit 

We pass now to consider the zero-temperature limit of 
the previous finite-temperature formalism. When T — > 0, 
the spacing 27rT between two consecutive Matsubara fre- 
quencies tends to zero. The sums over discrete Matsub- 
ara frequencies can then be replaced by integrals over 
continuous frequencies, provided the corresponding in- 
tegrands are not too singular (36[. We have then the 
equations 

Ep(fc) ^ J§G%(P-k)T(P), (12) 



r(p) = - 



II) r 



P 2 

2M 



J F (P) 



2ira y/2-K 

" (15) 

where we have defined M = m B + mp and /i = (ub + 
/xf)/2, while If(P), which results from the integration of 
the term with the Fermi function in Eq. ([5]) , is given by 
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(fc MF + k n y 



(*W + k P ) 2 - k 2 n 



kp 



kp 



(fc MF - k n y 



(10) where fc MF = v / 2tofMF; = tt-P while 



fen = (2m r )2 



P 2 

1- 2a + iQ. 

2M f 



i7rsgn(f2) 

(16) 

(17) 



Equations (llll) - (|16l) determine the thermodynamic 
properties of a Bose-Fermi mixture at zero temperature 
in the absence of boson condensation. They are thus rel- 
evant for a sufficiently strong coupling g, such that the 
system remains in the normal phase even at zero tem- 
perature. Upon lowering the coupling constant g, the 
condensation will start at a critical coupling g c , when 
the condition (fT0|) is first satisfied. 

We finally wish to remark the difference between the 
set of equations (TrTj) - ([T6|) . obtained by taking the zero- 
temperature limit of the finite-temperature formalism, 
and the equations that one would obtain, with the 
same choice of self-energy, by using directly the zero- 
temperature formalism (with the boson condensate set to 
zero) . In the latter case, one works on the real frequency 
axis and, more importantly, does not use the chemical 
potential up as an independent parameter, whose value 
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is fixed only at the end of the calculations by using the 
particle number equation. The chemical potential ap- 
pearing in the bare fermion Green's function remains in 
this way pinned to the noninteracting value. This makes 
the zero-temperature formalism less suitable when one 
wants to work beyond the purely perturbative regime, as 
in our case. 



III. FINITE-TEMPERATURE RESULTS 

The theoretical approach developed in section [TlJ can 
be used to explore the normal phase of a homogeneous 
Bose-Fermi mixture at arbitrary values of the boson and 
fermion masses and densities. In this section, where we 
present finite-temperature results, we focus on the spe- 
cific mass ratio me/mF = 87/40, relevant for the 87 Rb 
— 40 K mixture. More general mass ratios will be consid- 
ered at zero temperature. 
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2: (Color online) Critical temperature (in units of 

2 



FIG 

Ef = fc|/2mF) for condensation of bosons as a function of 
the boson-fermion coupling (fcpa) -1 for different values of the 
density imbalance (np — na ) / (np + iib) in a mixture with 
m B /ra F = 87/40. 



A. Critical temperature 

Figure [H presents the dependence of the condensa- 
tion critical temperature on the boson-fermion coupling 
(fcpa) -1 for a 87 Rb — 40 K mixture, at different values 
of the density imbalance (n-p — ^b)/("-f + «b)- The 
critical temperature was obtained by solving numeri- 
cally Eqs. ([3]) -([9]), supplemented by the condition ([TO]) , 
while the ending point at T = was calculated indepen- 
dently by solving the equations (fTTjl- (fl"6| (with the con- 
dition pop defining now the critical coupling g c at zero 
temperature) . The matching between finite-temperature 
and zero-temperature results confirms the validity of the 
equations derived in the zero-temperature limit, while 
providing simultaneously a check of the numerical calcu- 
lations. 

The overall behavior of the critical temperature as a 
function of coupling is similar to what found for equal 
masses 31|. The critical temperature starts from the 
noninteracting value (To = 3.31nj^ 3 /me) in the weak- 
coupling limit and eventually vanishes at a critical cou- 
pling g c , when the effect of the boson-fermion coupling 
is so strong that the Bose-Einstein condensate is com- 
pletely depleted, in favor of the formation of molecules. 
The weak dependence of the critical coupling on the den- 
sity imbalance previously found for equal masses is con- 
firmed also for this case with different masses. In this 
case, the critical coupling varies in the range 1.3 -j- 1.4, to 
be compared with the range 1.6-9-1.8 found for me = mp. 

A minor difference with respect to the case with equal 
masses is finally the presence of a weak maximum in the 
critical temperature, which reaches a value slightly above 
the noninteracting value To before its final decrease. This 
feature is more pronounced for intermediate values of the 
density imbalance. We attribute this feature to the del- 
icate balance between the effects of the boson-fermion 
interaction on the boson dispersion, and the condensate 
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FIG. 3: (Color online) Bosonic chemical potential at the crit- 
ical temperature T c as a function of the boson-fermion cou- 
pling (kpa)' 1 for different values of the density imbalance 
(rtF — tjb)/(^f +«b) in a mixture with 7TIb/?tif = 87/40. The 
corresponding fermionic chemical potential is reported in the 
inset. 



depletion due to molecular correlations. The boson dis- 
persion may in fact be hardened by the interaction (sim- 
ilarly to what one finds for a dilute repulsive Bose gas 
37| ) , thus leading to an increase of the critical tempera- 
ture. The predominance of one effect over the other one 
depends on a fine tuning of the mass ratio and density 
imbalance, and may explain the different behavior found 
for different values of these parameters. 



B. Chemical potentials 

Figure [31 reports the coupling dependence of the chem- 
ical potentials and /iF at the critical temperature for 
a 87 Rb — 40 K mixture at the same values of the density 
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imbalance considered in Fig. [2j 

In this case, the critical boson and fermion chemical 
potentials present the same qualitative behavior already 
found for a mixture with equal masses [3l|. The bo- 
son chemical potential decreases markedly with increas- 
ing coupling, and changes from /xb ~ 2irnpa/m r for weak 
coupling to fiB ~ — £o for strong coupling, with a small 
dependence on the density imbalance. 

The fermion chemical potential (reported in the inset) 
remains instead almost constant across the whole reso- 
nance. The decrease of the chemical potential due to the 
attractive interaction with the bosons is, in fact, partially 
compensated by the decreasing of the temperature when 
moving along the critical line (which increases fip) and 
by the Pauli repulsion between unpaired fermions and 
Bose-Fermi pairs. 



IV. ZERO-TEMPERATURE RESULTS 

The behavior of the critical temperature as a func- 
tion of the boson-fermion coupling discussed in the pre- 
vious Section, evidenced the presence of a quantum phase 
transition at zero temperature associated with a transi- 
tion between a supcrfluid phase with a boson conden- 
sate to a normal phase, where the condensate is com- 
pletely depleted. In this Section we examine in more de- 
tail this quantum phase transition by solving numerically 
the equations formulated at exactly zero temperature. 



A. Critical couplings and chemical potentials 
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FIG. 4: (Color online) Critical coupling g c as a function of 
the mass ratio (uib/wf) for different values of the density 
imbalance (np— tib)/(tif+^b). Panels (a) and (b) correspond 
to the ranges tob/wif < 1 and m^jm-p > 1, respectively. 



Figure [4] reports the critical coupling g c as a func- 
tion of the mass ratio tob/tof, with two distinct pan- 
els for the cases m^/m-F < 1 and 777.b/tof > 1. The 
different curves reported in Fig. 0] correspond to dif- 
ferent values of the density imbalance, ranging from 
the density-balanced case to the fully imbalanced one 
(np — ne)/(riF + tib) = 1-0, which represents the system 
with just one boson immersed in a Fermi sea. This is 
actually the same as a spin-down fermion sorrounded by 
a Fermi sea of spin-up fermions, since for a single particle 
the statistics is irrelevant. The critical coupling reduces 
thus to that for the polaron-to-molecule transition, re- 
cently studied in the context of strongly imbalanced two- 
component Fermi gases [3?j-|TH . The equations governing 
the one boson limit are reported in the Appendix. 

The critical coupling is strongly influenced by the mass 
ratio, especially for wb/wf < 1. In this case, for all 
values of the density imbalance, g c increases very rapidily 
as /tib/tof is decreased. For the single boson problem, 
an asymptotic expansion of the equations determining g c 
for tob/wf —> yields the result 

4 • 2 1 / 3 m F 32 • 2 1 / 3 , . 

9c « -= ~ + ( 18 ) 

07T TUB 107T 



« 0.535— +0.856 (19) 
m B 

which proves quite accurate even before the asympotic 
regime is reached (the deviation from the numerical so- 
lution is 15% for tob/tof = 1 and 2% for tyib /ttij? = 0.1). 
[See the Appendix for the details of the derivation of 
Eq. (HU).] 

The rapid increase of the critical coupling for the 
polaron-to-molecule transition when — !• was already 
noticed in Ref. l43l even though no asymptotic expression 
was reported there. Note however that, according to the 
analysis of the polaron-to-molecule transition of Ref. (4|| , 
for a mass ratio m^/m-j- < 0.15 the molecular state ac- 
quires a finite momentum in its ground state. Similar 
results were obtained in Ref. [4(| . The equations for the 
single boson problem here adopted assume that the for- 
mation of the molecule occurs at zero center of mass mo- 
mentum (as in Refs. 13814441) . The curve corresponding to 
the single boson problem in Fig. 3] may thus change for 
ms/ruF ^5 0.15, after taking into account the possibility 
of pairing at finite momentum. We expect however this 
change to be minor on the basis of our calculations with a 
finite boson density, which allow for pairing at finite mo- 
mentum and yield results close to the single-boson curve 
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also for me/niF < 0.15. 

Note further that the study of the three body system 
with two equal fermions with mass mp and one differ- 
ent particle with mass m B , shows that for m^/mp < 
0.0735 = (13.607) -1 the system is unstable due to a 
sequence of three-body bound states with energy — > 
— oo.\at\ A similar instability is expected to occur also 
in the many-body system with N equal fermions plus 
one different particle (a recent work has proven indeed 
that the above critical value for the three-body system 
provides a lower bound for the location of the instability 
in the many-body system [Hj]). This instability involv- 
ing three-particle correlations is out of the reach of our 
approach, which takes into account pairing correlations 
only. Also for this reason we have decided to limit our cal- 
culations to me/mF > 0.1 (for equal densities we stopped 
the calculations at me/mF = 0.2 because of numerical 
difficulties at small mass ratios due to the large variations 
of the bosonic chemical potentials for small variations of 
the masses and/or coupling values, see Fig. [5] below). In 
the regime where the three-particle instability is relevant, 
our equations, as well as the asymptotic expression (TT51) . 
will describe then a metastable situation, before the even- 
tual collapse of the system into a cascade of three-body 
bound states of increasingly lower energies. 

Consistently with our previous results, we observe a 
weak dependence of g c on the density imbalance. Such 
a weak dependence on the densities remains valid also 
for me > tof- In this case, all curves reach a minimum 
value of g c (=1.2-9-1.3) for mass ratios m B /mp in the 
range 3.5-^5, after which they increase slowly with the 
mass ratio. For the single boson problem with a large 
mass we have obtained the asymptotic expression (see 
Appendix) : 



g c « A(m B /m F ) 
where 



2 2/3 



1G 



1 



A(mB/mp) ' 37r A(tob/™f) 2 ^ ^ 



, . 2 4 / 3 / 4m B „ 
A(mB/m F ) = log 2 



(21) 



One can see from Eqs. (|2T))) and (|2"T|) that at large mass 
ratios g c increases very slowly, with a logarithmic depen- 
dence on the mass ratios. Due to this log-dependence the 
leading behavior, g c w (2 4 / 3 /7r) log ^f-, is reached only 
at extremely large mass ratios. Equation (|2"0")) , which in- 
cludes the first two corrections to the leading behavior, 
provides a better approximation, the deviation from the 
numerical solution being 15% for me/mF = 20 and 1.5% 
for me/mp = 100. 

Note finally that the two asymptotic expressions (p~8|) 
and (|21j) taken together imply the existence of a mini- 
mum value of g c at intermediate mass ratios, consistently 
with the results of Fig. 2J Interestingly, such a behavior 
of the single-boson problem captures also the overall be- 
havior of our curves at finite boson densities. 

At equal densities, our results for g c as a function of the 
mass ratio agree well with the results reported in Ref. 32 1 
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FIG. 5: (Color online) Critical coupling g c as a function of the 
density imbalance (hf — n^)/[rw + tib), for different values 
of the mass ratio ?tib / thjp . 



for the same case. This is because the equations used in 
Ref. [32| for calculating g c correspond, in a diagrammatic 
formalism, to the same choice of the self-energy as ours, 
but with the Dyson's equation expanded: G = Go + 
G0SG0, instead of G^ 1 = G^ 1 — E. Even though such an 
expansion is justified only when E is small , apparently 
it leads only to minor differences in the values for g c . 
For instance, for equal masses and densities we obtained 
g c = 1.62, to be compared with g c = 1.66 in Ref. 

According to the analysis of Ref. [32j, however, 



32|. 
only 

for sufficiently large values of the boson-boson scatter- 
ing length ae, the critical coupling g c lies in a stable 
region of the phase diagram. For small values of a B , the 
second-order quantum phase transition between a con- 
densed phase and a normal phase is in fact superseded 
by a phase separation between the two phases. According 
to our calculations, the compressibility matrix dn s /d(j, s > 
is always positive in the normal phase, indicating that 
the second order phase transition here explored lies, at 
worst, in a metastable region of the phase diagram. In or- 
der to examine its absolute stability within our approach, 
one should extend our study to the superfluid phase and 
make a comparison of the free energies for the normal and 
superfluid phases. This nontrivial extension is postponed 
to future work. We observe however, that if ratios a B /a 
of the order of 0.2-0.3 are sufficient to suppress phase 
separation in most of the phase diagram (as the results 
of Ref. [13] seem to indicate) , then the effect of a B on g c 
will be minor. 

The weak dependence of the critical coupling on the 
density imbalance is emphasized in Fig. [5j which presents 
g c as a function of the density imbalance for some rep- 
resentative values of the mass ratio tub/ nip. All curves 
show a weak dependence on the density imbalance, with 
a weak maximum at an intermediate imbalance (=0.55 
for equal masses, and similarly for the other mass ratios 
considered here). 

Figure[S]reports the chemical potentials at g c as a func- 
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FIG. 6: (Color online) Bosonic chemical potential at the crit- 
ical coupling g c as a function of the mass ratio ms/mF for 
different values of the density imbalance (up - tib ) / (wf + rie ) . 
The corresponding fermionic chemical potential is reported in 
the inset. 
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FIG. 7: (Color online) Bosonic momentum distribution curves 
at g c for a fixed density imbalance (wf — tib)/ (riF + riB) = 0.75 
and different values of the mass ratio (mB/mp). 



tion of the mass ratio tub/tof, for different density imbal- 
ances. As it can be seen from the main panel, the (nega- 
tive) bosonic chemical potential increases very rapidly in 
absolute value for ttib/tof — > 0. This is because at g c the 
bosonic chemical potential is already close to its strong- 
coupling limit, /^b ~ — £0) such that the dimensionless 
ratio \/1b\/Ef ~ 2mp / '{m r g 2 ) . Since g c ~ mp/mB for 
vtlb — > 0, we have \hb\/Ef ~ (mp/mB) 3 in this limit. 

In the opposite limit of large tub / mp , the ratio mp / m r 
slowly increases and eventually saturates to 1 for large 
mg, such that |/^b|/-Ef follows the slow logarithmic in- 
crease of g c in this limit. 

The fermionic chemical potential (reported in the in- 
set) depends weakly on the mass imbalance, reflecting 
the weak dependence on g c . As a matter of fact, the 
fermion chemical potential is determined essentially by 
the fermion density np , independently from the coupling 
value or mass ratio. 



B. Momentum distribution functions 

We pass now to study the momentum distribu- 
tion functions riB(|q|) and np(|k|) for the bosons and 
fermions, as obtained from Eqs. (fT"3|) and (fT4|) . respec- 
tively. We present results for a density imbalance (np — 
^b)/("f + "b) = 0.75, as to emphasize an interesting be- 
havior of the bosonic momentum distribution function, 
which occurs only for a sufficiently large imbalance. 

One can see, indeed, from Fig. [7] that the bosonic mo- 
mentum distribution function vanishes identically at low 
momenta. This empty region extends from |q| = up to 
a certain value |q| = go, which is determined essentially 
only by the density imbalance (for the specific case of 
Fig. Ego -0.55). 

The presence of the empty region can be interpreted 
as follows. For np 3> ub, most of fermions remain un- 



paired and fill a Fermi sphere of radius /ctjf — [i n F ~ 
^b)/67t 2 ] 1 / 3 , as Fig. |8] for the fermionic distribution 
clearly shows. At g c and larger couplings, the bosons 
are instead bound into molecules that, being compos- 
ite fermions, fill a Fermi sphere with a radius Pcf — 
(nB/671- 2 ) 1 / 3 . Now, as the region |k| < fcrjF is already 
occupied by the unpaired fermions, only fermions with 
|k| > fcr/F participate to the molecule. Since the momen- 
tum of the molecule is given by the sum P = q + k, the 
constraints |P| < Pcf and |k| > fcrjF then imply that 
only bosons with |q| > qq — fcuF — Pcf participate to the 
molecule, leaving thus empty the region |q| < qo- 

We have verified that the equation qo = ^uf — Pcf re- 
produces rather accurately the values of qo obtained nu- 
merically. In particular, the empty region does not exist 
at small density imbalance, when fcuF < Pcf- Note also 
that the initial rise of ns(|q|) after the threshold qo is due 
to the progressive increase of the phase-space volume cor- 
responding to the q's satisfying the above constraints at a 
given k. The saturation volume in phase-space is reached 
for |q| of the order of fcuF, after which ?^b(|<j|) starts 
to decrease, following essentially a molecular-like inter- 
nal wave-function riB(|q|) — |<^(|q| ) | 2 , where |</>(|q|)| 2 
can be approximated by the two-body normalized wave- 
function 0(|q|) = (87ra 3 ) 1 / 2 /(q 2 a 2 + 1), while the coeffi- 
cient bm can be interpreted as the molecular density. 

Even though these expressions are strictly valid only in 
the strong-coupling limit, when the binding energy eo S> 
-EuFi they work reasonably well already at g c . In particu- 
lar, the above expression for the molecular wave-function 
accounts for the differences in the curves calculated at dif- 
ferent mass ratios due to the strong dependence of g c , and 
then of a, on the mass ratio. Note further that the same 
molecular wave function describes the fermionic distribu- 
tion function at momenta |k| > fcuF, as it can be evinced 
from Fig. [5] (b). The comparison between Fig. [5] (b) and 
Fig. [7] shows indeed that the bosonic and fermionic dis- 
tribution functions become identical as the momentum 
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FIG. 8: (Color online) Fermionic momentum distribution 
curves at g c for a fixed density imbalance (np — riB)/(?iF + 
Ub) = 0.75 and different values of the mass ratio tyib / ttif ■ 
The two panels correspond to a different choice of the verti- 
cal scale. 



increases; at large momenta np(|k|) ~ 7iB(|k|) ~ C/k 4 , 
consistently with the universal large momenta behavior 
established in Ref. [49]. In particular, within our ap- 
proximation, one can prove by taking the large |k| limit 
in our expressions that the "contact" constant C is given 
by C = -Am 2 r f -£^zT(P)e iQ0+ . In the strong-coupling 
limit, where all bosons are bound into molecules, a com- 
parison with the expression for the molecular internal 
wave- function then leads to the equation C = Sttub/ol. 
Figure IHl reports the contact constant C at g c normalized 
to its strong-coupling limit value Stvub/cl, as a function 
of the mass ratio ms / mp for different values of the den- 
sity imbalance. One can see that the constant C at g c is 
close to its strong-coupling limit expression for all cases 
considered, with the largest deviations occuring at inter- 
mediate values of the mass ratio, as expected, since this 
is the region where g c reaches its minimum and conse- 
quently the strong-coupling condition is less respected. 
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FIG. 9: Contact constant C, normalized to its strong-coupling 
limit 8irnB/a, at g c vs. the mass ratio mB/mF for different 
values of the density imbalance (tif — wb)/(^f +%)• 



V. CONCLUDING REMARKS 

In this paper, we have studied a resonant Bose-Fermi 
mixture in the presence of a mass and/or density imbal- 
ance. We have first analyzed the finite temperature phase 
diagram for the specific case of a 87 Rb — 40 K mixture. 
We have found that the overall shape of the phase dia- 
gram is similar to what found previously for equal masses. 
The critical temperature starts from the noninteracting 

2/3 

value (To = 3.31n B /tub) in the weak-coupling limit and 
eventually vanishes at a critical coupling g Cl when the ef- 
fect of the boson-fermion coupling is so strong that the 
Bose-Einstein condensate is completely depleted. The 
critical temperature presents a weak maximum at inter- 
mediate coupling, not found for equal masses. We have 
explained this feature as resulting from the balance be- 
tween the hardening of the boson dispersion due to the 
boson-fermion interaction, and the condensate depletion 
due to molecular correlation. The predominance of one 
effect over the other one depends on a fine tuning of the 
physical parameters, thus explaining the different behav- 
ior found for different cases. 

We have then considered the zero-temperature limit. 
We have found that the critical coupling g c is signifi- 
cally affected by the mass imbalance, in particular for 
rne/rnF < f. In this case, for all values of the density 
imbalance, we have obtained a rapid increase of g c as 
tob/tof — > 0. We have also found that the density imbal- 
ance influences only weakly the critical coupling, as pre- 
viously found for equal masses [3l| . On the other hand, 
a sufficiently large density imbalance produces quite a 
remarkable effect on the bosonic momentum distribution 
function. We have found indeed that this function is 
completely depleted for momenta with magnitude be- 
low a certain value qo- This value is set by the differ- 
ence between the Fermi momenta associated with the 
unpaired fermions and the composite-fermions, respec- 
tively. Such a momentum distribution function reflects 
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the constraints on the internal molecular wave-function 
imposed by the presence of a large Fermi surface asso- 
ciated with the unpaired fermions. In this respect, the 
physics underlying this phenomenon is similar to that of 
the "Sarma phase" for a highly-polarized Fermi gas in the 
strong-coupling limit of the BCS-BEC crossover [50l - l55l |. 
In particular, the bosonic momentum distribution func- 
tion is similar to the momentum distribution function of 
the minority species in a polarized Fermi gas. Interest- 
ingly, the fact that in the present case the molecules can- 
not condense, but rather fill a Fermi sphere, smears some 
features of the momentum distribution function, but still 
preserves the presence of an empty region in the mo- 
mentum distribution function. Such an interesting effect 
could be measured in a Bose-Fermi mixture in the molec- 
ular limit by using the same technique used to measure 
the momentum distribution of an ultracold Fermi gas in 
the BCS-BEC crossover (H. 



nite boson density. The full knowledge of £r is however 
not necessary to calculate /ib in the limit ub — > 0. A 
study of the analytic structure of Eb in the complex fre- 
quency space along the lines of Ref. 43, shows that the 
limit riB — > corresponds, for g > g c , to the requirement 
that the minimum of the composite fermion dispersion 
fio(P) occurs exactly at zero frequency. By introduc- 
ing the retarded composite-fermion propagator via the 
replacement r R (P,w) = r(P,ifi ->• w + i0+), the dis- 
persion f2o(P) is determined by the pole of r R (P,cj) 
in the complex plane. By assuming the minimum of 
f2o(|P|) to occur at P = 0, one has then the equation 
r R (0, 0) _1 = 0, which determines hb at a given coupling 
and mass ratio. The critical coupling g c , on the other 
hand, is determined by the equation /ib = Eb(0, 0), as 
for finite density. In the limit ub — > only the pole of the 
fermionic Green's function contributes to the frequency 
integral in Eq. (fTTj) . yielding: 
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Appendix A: The one boson limit 



E B (0,0) = - 
where 

r R (P,£ F (P)) 



dP 



6R F (P)]r R (P,MP)), (Al) 
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1MB | 
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2m F M 



MP) 



(A2) 



In this Appendix we consider the "one boson" limit 
(rip — 7iB)/( n F + ^b) 1 of our equations. In this 
limit the fermions become free due to the vanishing 
boson density. One has then Ep = and fj, F = 
) 2 / 3 /(2m F ) = 2 2 / 3 E F . The boson self-energy E B 
remains instead finite and is determined by Eq. 
with no simplifications with respect to the case at fi- 



and 



MP) 



dp 



9KF(fP-p)] 

_P?_ _ 2iB £1 i |„ I 
2m r m F 2M ^ l^ B l 



A calculation of the above integral yields 



(A3) 
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4tt 2 
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where Rp = 
Lp = In 



pmrM-GPlmB/M) 2 ] 1 ^ while 
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2m r \nB\ 
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2tt 2 
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(A5) 

where 8m = m F — m B and m 2 — (m| + m B ) /2. By using 
the expressions (|A2I) to (|A5I) . the self-energy Eb(0, 0) can 



arctan 



few + |P| 

Rp 



arctan 



h — 22lF|p| 



Rp 



then be calculated easily by a simple one-dimensional 
integral over |P|. At g c the equation fiB = Eb(0, 0) then 
yields: 



Mb 



|^r R (|PU F (|P|)). (A6) 



The equation = r R (0,0)- 1 yields finally 
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2 \J m r |/i| arctan 
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(A7) 
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The simultaneous solution of Eqs. (|A6[) and (|A7[) allows 
to obtain g c and /is for a given mass ratio. Note finally 
that, even though we have derived Eqs. (|A6|) and (|A7j) as 
a limiting case of our equations for a Bose-Fermi mixture, 
they describe also the polaron-to-molecule transition in 
a Fermi-Fermi mixture, since in this limit the statistics 
of the minority species becomes immaterial. We have 
verified, indeed, that our results for q r in this limit agree 
with the results reported in Ref. |43j |. with the same T- 
matrix choice of the self-energy, for some specific mass 
ratios. 



Asymptotic expressions for c 
mass ratios 



at small and large 



We conclude this appendix by presenting the deriva- 
tion of the asymptotic expressions (fT8j) and (|20|) . We 
assume /ib to be large and negative as it occurs when g 
is large. The validity of this assumption is verified by 
the asymptotic expressions for g c that are obtained ac- 
cordingly. The chemical potential [i = (^b + Mf)/2 is 
then also large and negative. By expanding Eq. (IA7|) in 
powers of / \/m r \fj,\ one obtains 

4 



2fi « -e 



K a. 

3irm r 



(A8) 
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-k't a in Eq. 



Note that the subleading term 
describes the mean-field repulsion between the molecule 
forming in the strong-coupling limit and the fermions, 
as it can be seen by casting it in the form n F „f^ F «mF; 
where wmf = Mm-p/ (tof + M) is the reduced mass of a 
molecule and one fermion, while 



a MF 



(1 + m F /m B ) 2 
1 + tof/2tob 



(A10) 



is the molecule-fermion scattering length within the Born 
approximation [57|. 

The large and negative value of /zb then implies 
/f(P) ~ "-f/ImbI ~ "-f/eo: a s it can be seen more easily 
directly from Eq. (|A3|) . The expansion of (|A2|) for |^b| 
large then yields 



r R (P,? F (P)) 



where we have used Eq 
using Eq. (IA1[) we obtain 
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, and a = mB/(2Afmp). By 
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FIG. 10: (Color online) (a) The asymptotic expression (|A15|) 
for g c at small mass ratios (dashed-dotted line) is compared 
with the full numerical solution (full line) . The dotted curve is 
obtained by neglecting the subleading term within the brack- 
ets in (|A15[l . (b) The asymptotic expression (|A20|l for g c at 
large mass ratios (dashed-dotted line) is compared with the 
full numerical solution (full line). The dashed curve is ob- 
tained by neglecting the second and third terms on the right- 
hand side of (|A20j) . while the dotted curve neglects only the 
third term. 



By expanding the expression (|A12|) for a small mass ratio 
me/niF *C 1, we obtain 



S B (0,0) 



1 mg 

3 TTlF 

2 4 / 3 k F m F 



2_ / mg 
15 \ mF 
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37r aniB ttib 
The equation /ib = Sb(0, 0) then yields 
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5 tof 



2 4 / 3 k 
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(A13) 



(A14) 



where we have kept only the leading term in the expres- 
sion (| A9|) for [Ib (the term /ip would give a correction 
of order (tob/tof) 2 to the asymptotic expression for <? c , 
with the last term in (|A9|) contributing an even smaller 
correction). By substituting eo = l/(2TO r a 2 ) we obtain 
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then 



1 



4 • 2 1 / 3 m F 
37r m B 



8 m B 

5 TOp 



(A15) 



which coincides with the expression (|18p reported in 
Sec. IV, and gives a large value of g c for small mass 
ratios, consistently with our starting assumption. The 
asymptotic expression (|A15[) is compared with the full 
numerical calculation of g c in Fig. [10] (a) . 

The expansion of Eq. (|A12[) for a large mass ratio 
m B /mp 3> 1 yields instead 



E B (0,0) 
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2 4 / 3 fc F 
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1 m F \ 
2 l0g 4^J 
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where we have disregarded corrections to (|A17[) smaller 
at least by a factor mF / ro B ■ By equating (|A17|) to (|A9[) 
we get then 
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from which, by multiplying both sides of Eq. (|A18[) by 
2amp/fc F , we obtain 

1 21/ 4m B \ 16(fc F a) 2 , A „. 

— = — log — - - 2 -23fc F a+ V (A19) 
kpa 7r \ mF / dw 



which shows that, to leading order, g c w log ^ at 
large mass ratios. The slow log-dependence of g c on 
m B / mF makes however significant to keep also a few sub- 
leading corrections to the leading behavior. The solution 
of Eq. (|A19|) by iteration yields then 



kpa 7r V mp / 2 



25 



4m b 



- 2 
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16 
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3tt 



4 (log 



4m B _ 2 



(A20) 



which ignores corrections of order l/(log ^-) 3 or higher, 
and coincides with the expression (f2"U|) reported in 
Sec. IV. The asymptotic expression (|A20[) is compared 
with the full numerical calculation of g c in Fig. [TO] (b). 

As a final remark, we note that it is the increase of the 
self-energy S B (0, 0) both at small and large mass ratios 
to drive the increase of g c in the two limits. Since the po- 
laron energy obtained with our self-energy coincides with 
that obtained with the variational wave-function intro- 
duced by Chevy in [5811 . which is deemed quite sensible 
for the polaron energy |43f , we believe our results for the 
leading behavior of g c (and thus for the polaron-molecule 
transition) at small and large mass ratios to be quite ro- 
bust. 
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